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q ■ Abstract 

We study quantized beamforming in wireless amplify-and-forward relay-interference networks with 
any number of transmitters, relays, and receivers. We design the quantizer of the channel state infor- 
mation to minimize the probability that at least one receiver incorrectly decodes its desired symbol(s). 
Correspondingly, we introduce a generalized diversity measure that encapsulates the conventional one 
as the first-order diversity. Additionally, it incorporates the second-order diversity, which is concerned 
with the transmitter power dependent logarithmic terms that appear in the error rate expression. First, 
we show that, regardless of the quantizer and the amount of feedback that is used, the relay-interference 
network suffers a second-order diversity loss compared to interference-free networks. Then, two different 
quantization schemes are studied: First, using a global quantizer, we show that a simple relay selection 
scheme can achieve maximal diversity. Then, using the localization method, we construct both fixed- 
length and variable-length local (distributed) quantizers (fLQs and vLQs). Our fLQs achieve maximal 
first-order diversity, whereas our vLQs achieve maximal diversity. Moreover, we show that all the 
promised diversity and array gains can be obtained with arbitrarily low feedback rates when the 
£ i transmitter powers are sufficiently large. Finally, we confirm our analytical findings through simulations. 
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Index Terms 

Wireless relay network, beamforming, interference, distributed vector quantization, symbol error 
probability, diversity gain, array gain. 

I. Introduction 

While it has been demonstrated in several studies that cooperation can greatly improve the 
performance and reliability of wireless network communications [HJ — [O, interference still re- 
mains to be a fundamental issue in cooperative network design. Most of the previous work on 
cooperative networks relies on orthogonal channel allocation so that different transmitters do 
not interfere with each other. However, allocating orthogonal channels for each user may not 
be desirable due to time and bandwidth limitations |[6l, 0. In such cases, one should explore 
effective ways to deal with interference while preserving cooperative diversity gains. 
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Multiple antenna interference cancelation techniques are very effective when dealing with 
interference in cooperative networks [8J. They offer reasonable performance with low decoding 
complexity. In this work, we consider a different approach. To be able to study the ultimate 
performance limits, we do not put any restrictions on our decoders. We would like to design a 
cooperation scheme that achieves maximal diversity benefits, and thus provides high reliability, 
even in the presence of multiuser interference. 

For networks with a single transmitter-receiver pair and no interference, network beamforming 
using amplify-and-forward (AF) relays has shown to achieve the maximal spatial diversity 
BH, ifTOl . However, the optimal beamforming policy requires one or two real numbers to be 
broadcasted from the receiver to the relays. Using distributed beamforming with quantized 
instantaneous channel state information (CSI), it is possible to obtain both maximal diversity, 
as well as high array gain with only a few feedback bits from the receiver [fTTTl - [fT3l . A special 
case of quantized feedback for cooperative networks is the relay selection scheme |[T4 | - lfT6ll . It 
has been formally shown in |fTTf that, for a network with R parallel relays, the relay selection 
scheme provides the maximum diversity R. 

Quantized feedback schemes have also been studied for non-cooperative multiuser interference 
networks. In iTPTll . the author considers zero-forcing beamforming with finite rate feedback in 
multiple-input multiple-output (MIMO) broadcast channels. Interference alignment for multiuser 
interference networks with limited feedback has been studied in [18]. Unlike what we shall study 
in this work, where we seek to optimize the reliability of the system in terms of the diversity 
gain, the goal of the above two papers was to optimize the data transmission rate in terms of the 
multiplexing gain. A common conclusion that we can infer from both studies is that, in order to 
achieve the same multiplexing gain as a system with perfect CSI, the feedback rate should be 
increased at least logarithmically with the transmitter power; any constant feedback rate results 
in a complete loss of multiplexing gain. This is unlike point-to-point systems where feedback 
is not even necessary to achieve the maximal multiplexing gain [TTvTl . and a few feedback bits 
is usually sufficient to transmit with rates that are close to the one with perfect CSI lU9ll . The 
feedback requirements of interference networks appears to be considerably higher than that of 
interference-free networks. 

What are the feedback requirements if instead we would like to ensure maximal reliability in 
the presence of interference? One goal of this paper is to answer this question for cooperative 
networks with K transmitters, L receivers, and R parallel AF relays. We assume that each 
transmitter and each relay has its own short term power constraint. The transmitters do not have 
any CSI. Each receiver knows its own receiving channels and the channels from the transmitters 
to the relays. Each relay only knows the magnitudes of its own receiving channels. Each relay 
and each receiver also has partial CSI provided by feedback. The feedback information represents 
a quantized beamforming vector. In that sense, this paper is also a generalization of single-user 
quantized network beamforming ifTTTl to multiuser interference networks. On the other hand, 
such a generalization is quite challenging because of the distributed nature of the network. Let 



us now describe some of these challenges and our approaches to address them. 

In interference networks, the relays amplify both noise and interference, which results in 
completely different problem formulations and solutions. Second, there are multiple receivers 
that have different optimal beam directions. As a result, it is difficult to design a scheme that 
can provide a reasonable performance to all the users. 

Another difficulty is related to acquiring feedback information from several separated receivers. 
The optimal beamforming policy requires the full CSI of the interference network. In practice 
however, none of the receivers can obtain such information via training methods. We thus 
consider two different quantization schemes: In the first scheme, the feedback information is 
provided by a global quantizer (GQ) that knows the entire CSI. We use this hypothetical quantizer 
to analyze the performance limits of network beamforming in the presence of interference. In the 
more practical second scheme, we use distributed local quantizer (LQ) encoders at each receiver. 
Each receiver can access only a part of the CSI, and provides its own feedback information for 
the relays and the other receivers. 

In EOl . we introduced a general systematic LQ design method, called localization, in which 
one synthesizes an LQ out of an existing GQ using high-rate scalar quantization combined with 
entropy coding. In the same work, we described an application of the method to MIMO broadcast 
channels. In this work, we apply it to design LQs for our network model. Therefore, our GQ 
has another important purpose other than the one we have previously mentioned: It will also 
serve as the basis of our LQs. 

We would also like to note that the LQ design in this paper distinguishes itself from the one 
in EOl in several ways, even though the underlying localization method will be the same. First, 
we need to consider a totally different and much more complicated distortion function. Second, 
the high-rate scalar quantizers, that form the crucial part of the method, should be designed 
accordingly. Third, the performance analysis of the resulting LQs is thus different and more 
complicated. As a result, in this work, we will only analyze the performance of localization for 
a particular class of GQs that are based on relay selection. 

Our performance measure is what we call the network error rate (NER). Given a fixed channel 
state, it is the probability that at least one user incorrectly decodes its desired symbol(s). In that 
sense, any receiver can be interested in the symbols transmitted by any subset of transmitters. 

We use a generalized diversity measure to characterize the asymptotic behavior of the NER as 
the transmitter powers grow to infinity. In what follows, we describe this measure together 
with its motivations: Suppose that a wireless communication system achieves an error rate 
of C(P a log' 3 P) _1 , where P is the transmitter power constraint and C is a constant that 
is independent of P. Then, we call a and f3, the first-order and the second-order diversity 
gains, respectively, and say that the scheme achieves diversity (a, (3). Such a definition of 
diversity is more precise than the traditional one as we demonstrate by an example: For two 
hypothetical communication systems with diversity gains (a, fix), and (a,/3 2 ), where a > 1 and 
oo > Px > f3 2 > — oo, the former always outperforms the latter for all P sufficiently large. On 



the other hand, the traditional definition, according to which the diversity gain is a for both 
systems, fails to distinguish between the asymptotic performance of the two. 

The main contributions of this paper can be summarized as follows: First, we show that, 
regardless of the quantizer and the amount of feedback that is used, the maximal achievable 
diversity of our network model is (R, —R) when K > 1, whereas it is (R, 0) when K = 1 U In 
other words, the relay-interference network suffers from a second-order diversity loss compared 
to an interference-free network that can achieve diversity (R, 0) with K — L — 1 [fTT|. Then, 
we construct a relay-selection based fixed-length GQ (fGQ) that can achieve maximal diversity 
for any K. Next, using our fGQ and the localization method, we design both fixed-length and 
variable-length LQs (fLQs and vLQs). Our fLQs can achieve diversity (R, —2R) when K > 1, 
and diversity (R, —R) when K = 1, using R feedback bits per receiver. They show that it is 
possible to achieve very high reliability using a fixed number of feedback bits. On the other 
hand, our vLQs can achieve maximal diversity gain for any K. Moreover, the feedback rate they 
require decays to zero as the transmitter powers grow to infinity. Therefore, they provide a very 
fortunate answer to the question that we have posed earlier: In a relay-interference network, it 
is possible to achieve maximal reliability using arbitrarily low feedback rates per receiver, when 
the transmitter powers are sufficiently large. Another desirable property of our vLQs is the fact 
that the array gain they provide can be made arbitrarily close to the one provided by the fGQ. 

The rest of the paper is organized as follows: In Section HH we introduce our network model, 
performance and diversity measures, and problem definition. In Section [Till we show that the 
maximal diversity of our network model is (R, —R). In Sections [IV] and |Vj we introduce our GQ 
and LQ designs, respectively. Numerical results are provided in Section |VT1 In Section IVIH we 
draw our major conclusions. An upper bound on the probability density function (PDF) and the 
cumulative distribution function (CDF) of a frequently used random variable (RV) is provided 
in Appendix [A] Some other technical proofs are provided in Appendices |B] through E 

Notation: For a logical statement S, "S is true for x sufficiently large" means that there exists 
x < oo such that for all x > x , S is true. || • || indicates the 2-norm, || • H^ is the infinite 
norm, (-|-) is the inner product. C, K. and Z + represent the sets of complex numbers, real 
numbers, and positive integers, respectively. det(A) is the determinant of a square matrix A. 
A T , A H denote the transpose and the Hermitian transpose of A, respectively. P represents the 
probability. f x (-) is the PDF, and F x (-) is the CDF of an RV X. E[X] is the expected value 
of X. X ~ T(k,9) means that X is a Gamma RV with fx{%) = x ^w^ f° r x > and 
fx{x) = for x < 0, k, 9 > 0. For any sets A and B, A — B is the set of elements in A, but not 
in B. \A\ is the cardinality of A. A r = {(ai, . . . , a r ) : a 1; . . . , a r G A}, r E Z + , is the cartesian 
power. 7 e = 0.577... is the Euler-Mascheroni constant, e = exp(l), and is the empty set. For 
a real-valued function / : C -)• R with C C C K , let M = {x : x e C, /(x) = max x / eC /(x')}. 

'The case K = 1 corresponds to a relay-broadcast network that does not suffer any multiuser interference. Even though our 
main goal in this paper is to analyze interference networks, we present the extension of our results to broadcast networks, so as 
to demonstrate the detrimental effects of interference in a comparative manner. 



Then, arg max xeC /(x) is the unique vector x* with the property that x* -< x, Vx 6 .M, and "-<" 
represents some partial ordering (e.g. lexicographical ordering) of complex vectors. We define 
argmin(-) in a similar manner. Finally, log(-) is the natural logarithm, log 2 (-) is the logarithm 
to base 2, cosh(-) is the hyperbolic cosine, Q(-) is the Gaussian tail function, T(-) is the gamma 
function, E±(x) = f™ e~ 1 e~ xt dt is the exponential integral, and K u (-) is the modified Bessel 
function of the second kind of order v. 

II. Network Model and Problem Statement 
A. System Model 

The block diagram of the system is shown in Fig. [Q We have a relay network with K 
transmitters, L receivers, and R parallel relays. The cases K — 1 and K > 1 correspond to a 
relay-broadcast network and a relay-interference network, respectively. We assume that there is 
no direct link between the transmitters and the receivers. 
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Fig. 1: System block diagram. In the figure, TX fc , Relay r , and RX^ stand for the A;th transmitter, 
rth relay, and the £th receiver, respectively. 



Denote the channel from the A;th transmitter to the rth relay by fk r and the channel from 
the rth relay to the fth receiver by g r £. Let h = (/ n , . . . , f KR , #n, . . . , g RL ) denote the channel 
state of the entire network. We assume that the entries of h are independent and distributed as 
fkr ~ £/V(0, cr 2 f kr ), 9re ~ CM(0, a 2 gri ) with finite variances cr fkr , a 9rl < 00, Vr, k, L For brevity, 
let g£ — (gu, . . . , giu), which denotes all the channels from the relays to the £th receiver. 

Only the short-term power constraint is considered, which means that for every symbol 
transmission, the average power levels used at the fcth transmitter and the rth relay are no 
larger than P Sk and Pr t , respectively. 

We assume a quasi-static channel model; the channel realizations vary independently from 
one channel state to another, while within each channel state the channels remain constant. We 



assume that the £th receiver knows gi and each relay knows the magnitudes of its own receiving 
channels, i.e. the rth relay knows \fk r \, k = 1, . . . , K. Some possible procedures to reveal the 
channel states to the receivers can be found in ifTTTl . lfT3l . For completeness, we give an outline 
of one possible way: The £th destination can acquire the knowledge of g r e by training from the 
rth relay. The rth relay can acquire the knowledge of \fk r \ using training sequences from the 
fcth source. It can also amplify and forward its received training signal from the source to the 
destination, so that the destination can estimate the product of fk r and g r £. As g r £ is known by 
the destination, /& r can be estimated. 

Each relay and each receiver also has partial CSI provided by feedback. In this paper, we 
consider two different feedback schemes, namely the global and local quantization schemes. 

B. Global Quantization 

Our global quantizer GQ is defined by a global encoder and a global decoder, as described 
in Fig. [2l The global encoder consists of two parts. For each channel state, first, a GQ encoder 
QGE : C R(X+L) -> X G maps the channel realization h to an index in X G = {1, . . . , |X G |}, the 
index set of the codebook elements. Then, a lossless global compressor GQC : X s — > J G maps 
this index to a binary description. 
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Fig. 2: Global quantizer operation. 



Let l(j) denote the length of a binary description j. We call GQ a fixed-length GQ (fGQ) if 
l(j) = [log 2 |X G |], Vj E J G . Otherwise, we call GQ a variable-length GQ (vGQ). 

In either case, the global encoder feeds back GQC(GQE(h)), using l(GQC(GQE(h))) bits. The 
feedback bits are received by the global decoders without any errors or delays. 

There is a unique global decoder at each relay and each receiver, which comprises of the com- 
plementary parts to the global encoder: A lossless decompressor and a quantizer decoder. First 
the decompressor GQC -1 : J Q — Y X G reconstructs the quantization index from the received binary 
description. It is followed by the quantizer decoder GQD : X G — > C G which maps the quantization 
index to a codebook element. The codebook C G has |X G | elements, C G = {xi, . . . , xixgi}. Without 
loss of generality, for GQE(h) = m, we set GQD(m) = x m e C 9 . For the rest of this paper, we 
will use the well-known notation Q 9 (h) = (GQD o GQC -1 o GQC o GQE)(h) = (GQD o GQE)(h). 
Therefore, GQ : C R{K+L) -> C G , and GQ(h) = x, for some x e C G . 

In the most general case, the rth relay may make use of the side information \f kr \ in the 
process of decoding the feedback information. However, in order to keep the relay operation as 
simple as possible, we do not consider such a scenario in this paper. 



C. Local Quantization 

We define our local quantizer LQ by L local encoders, with the ith encoder at the ith receiver, 
and a unique local decoder at each receiver and relay, as described in Fig. [3] The £th local 
encoder comprises of two parts: An LQ encoder LQE^ : C R ( A+1 ) — > X\ and a lossless local 
compressor LQC^ : X\ — > J\. Note that the domain of each LQ encoder is different from the 
domain of the GQ encoder. For the £th encoder, the domain corresponds to the channel states 
from the transmitters to the relays and from the relays to the £th receiver, represented by the 
concatenation vector [f , g^] . 
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Fig. 3: Local quantizer operation. 



The M receiver feeds back LQQ(LQE^([f, g^])), using [(LQQ(LQE^([f,g^]))) bits. We call LQ 
an fLQ if, l(j) = |~log 2 \Jl\~\, Vj G J\, W. Otherwise, we call it a vLQ. For the latter case, the 
feedback rate of the £th receiver can be expressed as R^(LQ) = E[[(LQQ(LQE^([f, &])))]. 

After all the L feedback messages are exchanged between the receivers and the relays, each 
of them decodes the feedback bits using the local decoder. The local decoder is the composition 
of a decompressor LQCT 1 : \\ t J\ — > Yle^e an d a quantizer decoder LQD : Yli^i ~* C L - Overall, 
LQ(h) 4 LQD(LQE 1 ([f,g 1 ]),...,LQE L ([f,g L ])). Thus, LQ : C^+ L ) -► C\ and LQ(h) = x, for 
some x G C L . 



D. Transmission Scheme 

We use a two-step AF protocol iflOl . [fTT]]. In the first step, the A;th transmitter selects a symbol 
Sk from a constellation 5%, where \Sk\ < oo, P(s fe ) = |5fe| _1 , Vs^ G <S&, and sends ^jPs k Sk- We 
normalize Sk as E[|sfc| 2 ] = 1. Thus, the average power used at the A;th transmitter is Ps k . During 
the first step, there is no reception at the receivers, but the rth relay receives 

K 

U = ^ fkrSk\fPs~ k + VOr, (1) 

fc=l 

where r] 0r ~ CN(0, 1). 

Suppose that a quantizer Q : C R ( K+L ) — > C, global or local, is employed in the network, 
and Q(h) = x, for some x G C. Then, the relays use the beamforming vector x to adjust their 



transmit power and transmit phase. During the second step, the transmitters remain silent, but 
the rth relay transmits 

Ujf JU f \f L/<y by* . \ / 

where the relay normalization factor p r is given by 

Pr = J? • (3) 

The average power used at the rth relay can be calculated to be E Sl: ,„ :SK>Vor [\u r \ 2 ] = \x r \ 2 PR r , Vh. 
We require < \x r \ < 1 as a result of the short term power constraint. The channel state 
dependent normalization factors p r ensure that the instantaneous transmit power of each relay 
remains within its power constraint with high probability!] 

Also, note that within the restriction of < |x r | < 1, p r is the maximal normalization factor 
that we can use. In other words, if a factor p" r satisfies p" r > p r for some h, then it violates 
the short term power constraint. Still, one can employ another factor p' r with p' r < p r , Vh (e.g. 
p' r = Pr t /(1 + X]fc=i(l + \fkr\ 4 )Ps k ))- We shall discuss later in Section HIT1 whether or not such 
a different choice of the normalization factor can improve the network performance. 

After the two steps of transmission that has been described above, the received signal at the 
£th receiver can be expressed as: 

K R R 

Ve. = ^ 5Z X ryfp~rfkr9re\fPs i Si + 22 X r9re^/~Pr~ri0r + VU, (4) 

k=l r=l r=l 



where rju ~ CN(0, 1) is the noise at the £th receiver. We assume that the noises r] 0r , r — 1, . . . , R, 
and r/u, £ — 1, • • • , L are independent. 

E. Performance Measure 

The £th receiver attempts to decode the symbols of the transmitters with indices given by an 
arbitrary but fixed set T>t C {1, . . . , K}, Vi ^ 0. As an example, for a network with K = 3 
and L = 2, let V>\ = {1,2} and T> 2 = {2,3}. Then, the first receiver is interested only in the 
symbols of the first and the second transmitters, while the second receiver is interested only in the 
symbols of the second and the third transmitters. In general, we assume that {J f U e = {l,..., K}. 
This guarantees that at least one receiver is interested in the symbols of the /cth transmitter. In 
particular, for K = 1, we have T>e = {1}, W. 

Let us call the vector of transmitted symbols s^ = [sk\k&T> e as the super-symbol relevant to the 
£th receiver, and s^ be its decoded version. We say that an error event occurs at a receiver if it 
incorrectly decodes its desired super-symbol. In this case, the optimal decoder at the £th receiver 

2 Because of the noise at its received signal, a relay can exceed its transmit power constraint at some instants. The phrase 
"short-term" comes from the observation that, regardless of the channel states, the relay always obeys its power constraint when 
its transmit power is averaged over the transmitted symbols and the noise. 



is an individual maximum likelihood (ML) decodeo given by s^ = arg max s ' e ^ P(s' e \y e , x, h), 
where S^i = Ylkev $k is the relevant super-symbol alphabet. For a fixed channel state h, and 
beamforming vector x, let SER™ L (x, h) = P(s£ ^ sg) denote the conditional super-symbol error 
rate (SER) of the £th receiver with the individual ML decoder. 

Let us now define a single quantity that represents the SER performance of all the receivers. We 
define the conditional network error rate (conditional NER, or CNER), denoted by CNER(x, h), 
as the probability that at least one receiver incorrectly decodes its desired super-symbol. 

Our performance measure, the NER, is the expected value of the CNER. Given a quantizer Q 
global or local, the NER can thus be expressed as 

NER(Q) = E h [CNER(Q(h), h)]. (5) 

F. Diversity Measure 

Let us also define a unique diversity measure for our network. Let P Rr = Pr,,P, r = 1, . . . , R, 
P Sk = Ps k P, k — 1, . . . ,K, where Ps k ,Pn, r < oo. In other words, we allow the power constraint 
of each transmitting terminal to grow linearly with P. Then, the first-order diversity achieved 
by a quantizer Q is given by 

dM A „,„ -M, (6) 

P^oo logP 

One problem with this conventional definition of diversity is that it fails to characterize the 
asymptotic effect of possible sub-linear P-dependent terms (e.g. logarithmic terms) in the error 
rate expression. In order to properly handle such cases, we define the second-order diversity as 

& lim j a g«a<q) + d,(q)iogP j 

P^oo log log P 

Note that the first-order diversity is always positive, while the second-order diversity can be 

negative. 

Now, the diversity (gain) achieved by a quantizer Q is given by d(Q) = (di(Q), ^(Q))- 
With these definitions, the asymptotic performance with a quantizer Q, as P grows to infinity, 

can be expressed as 

NER(Q) = G A (P)(logP)- d2(Q) P^ l(Q) , (8) 

where the factor G A (P) is the array gain. It is sublogarithmic in the sense that limp^oo j^p = 0. 
Also, we use it only when we compare the performance of two quantizers that provide the same 
diversity gain. 

Finally, for two diversity gains d — (di, d 2 ), and d! = (d[, d' 2 ), we say that d is higher than d' 
(or d > d') if either d\ > d[ or d 1 = d[, d 2 > d' 2 . 

3 In the literature, the phrase "individual" usually refers to the cases in which the a posteriori probability is maximized over 
a single transmitter alphabet. Note that, in our case, the maximization is over the product alphabet Se that represents the set of 
all super-symbols that the Ah receiver is interested in. 
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G. Problem Statement 

Our goal is to design the quantizer Q, given a limited feedback rate, such that the NER is 
minimized. We consider this problem for both GQs and LQs. 

To achieve our goal, we first determine the maximal possible diversity with our network model. 
Then, we design structured fGQs that can achieve this diversity. Finally, we use our observations 
on fGQs to systematically design fLQs that achieve maximal first order diversity, and then, vLQs 
that achieve maximal diversity. 

We would like to note that, as demonstrated in |[TTT|. the numerical optimization of our 
quantizers is always possible by using algorithms such as the Generalized Lloyd Algorithm 
fl2TJ, ll22ll . These algorithms can be used to improve the array gain performance, or in some 
particular cases, the second-order diversity performance of our structured codebook designs. We 
will not consider such optimizations in this paper since they are straightforward. 

III. Lower Bounds on Quantizer Performance 

Before we attempt to design a high-performance low-rate quantizer, it is natural to determine 
the best possible performance we can expect with any quantizer. In this section, we find lower 
bounds on the NER for both relay-interference and relay-broadcast networks that hold for any 
quantizer Q, global or local. 

Let A" = {x6 C^ : || x Hoc < 1} represent the set of all beamforming vectors. Then, we have 

Theorem 1. Let Q : C /?(x+L) -> C with C C X. Then, there are constants <C 1 ,C 2 < 00 that 
are independent of both P and Q, such that for all Q, and for all P sufficiently large, 

NER(Q) > d-^, K = l, 

\og R P (9) 

NER(Q) > C 2 -|^, K>\. 

Moreover, the bounds in ((9]) hold for any relay normalization factor p' r that satisfies p' r < p r , Vh. 

Proof: Please see Appendix E ■ 

In other words, for relay-broadcast networks, the maximal diversity gain is (R, 0). Indeed, for 

a network with K — L — 1, it was shown in [fTT]| that diversity (R, 0) is achievable. 

On the other hand, for relay-interference networks, the maximal diversity gain is (R,—R). 

Since (R, 0) > (R,—R), interference results in a second order diversity loss in our network 

model. 

Theorem [Q also shows that a different relay normalization factor p' r cannot improve the 

diversity upper bounds, provided that it satisfies the short-term power constraint, and a codebook 

C C X is employed. Thus, for the rest of this paper, we will only consider p r as our relay 

normalization factor. 
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An immediate question that stems from Theorem Q] is whether there exists finite rate quantizers 
that can achieve maximal diversity. In the next section, we construct an fGQ that provides an 
affirmative answer. 

IV. Maximal Diversity with an fGQ 

In order to determine an fGQ that can achieve maximal diversity, let us first determine, for 
any K, the optimal GQ given a fixed codebook with finite cardinality. 

Proposition 1. Given a fixed codebook C with \C\ < oo, the optimal GQ is given by GQ£(h) = 
argmin xe c CNER(x, h). 

Proof: Let Q' : C R ( K+L ^ -> C. We have 

CNER(GQc(h),h) < CNER(Q'(h),h) =*► NER(GQc) < NER(Q'), VQ'. (10) 

Thus, GQ£ performs at least as good as any quantizer with codebook C. ■ 

Therefore, given that we employ an optimal GQ encoder given by Proposition [Q, the GQ 
codebook uniquely determines the system performance. But, there is one complication: If we 
ever want to implement the optimal GQ encoder, we should be able to evaluate CNER(x, h), 
for any given x and h. Unfortunately, a closed form characterization of the CNER is very 
difficult, if not impossible. For that reason, we design a suboptimal quantizer that, instead of the 
actual CNER, uses an upper bound on the CNER. Fortunately, this suboptimal quantizer will be 
powerful enough to achieve maximal diversity for any K. 

A. An Upper Bound on the CNER 

For the £th receiver, instead of the individual ML decoder described in Section III-EL suppose 
that we employ a joint ML decoder s^ = argmax s / e y P(s'|^, x, h), where 5? = Y\ k S^- Recall 
that, for the individual ML decoder at the £th receiver, the a posteriori probability was maximized 
over Yliev •%• F° r the joint ML decoder, the maximization is over Ha:^ at a ^ tne receivers. 

Let SERf L (x, h) = P(s^ ^ s) denote the error rate of the joint ML decoder. Then, we have 
SER™ L (x,h) < SERf L (x,h), W. Also, from ©fl 

SERf L (x,h) <^ Q( v /2 7 £, s ,s(x,h)), (11) 

where 



k=1 {S k ~ S k )^F Sk }^ r= i Jkr^/Pr9HX r \ 

7£ S s(x, ni = ^ . (12) 

4 U +Ljr=lPr\9rer\Xr\) 

Note that, in order to be able to perform ML decoding, the receivers should know which beamforming vector is used by the 
relays. In other words, for each h, the receivers should know Q(h). This explains why we need to have a quantizer decoder at 
each receiver as well as each relay. 
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In (fTTI) and (fT2~l) . the decoded symbol vector for each receiver is obviously different, i.e. Sg, 
though we have omitted the dependence on £ for brevity. Furthermore, from now on, we shall 
omit the condition s, s G 5? in the summations as it is clear from the context. 

Now, using a union bound over all the receivers, it follows for the CNER that 

L 

CNER(x, h) < J2 SERf L (x, h) (13) 

L 

<]TsERf L (x,h) (14) 



! L 

' ' 1=1 s^s 



(15) 



This upper bound can easily be evaluated for any constellation and thus, it is good enough for 
our purposes. However, for clarity of exposition in the rest of the paper, we seek a much simpler 
bound. First, let us define 

7* (x, h) = min 7^ (x, h) , (16) 

and 

7 L (x, h) = min 7^(x, h) (17) 

= min min 7^ (x, h). (18) 



Then, (|T5l) can be further bounded as 

CNER(x,h) < ^^Q(V2 7f (x,h)) (19) 

£=1 

< L(l ^~ 1} max Q ( v/2 7 £(x,h)) (20) 

= 2C Q( v / 2 7 L (x,h)) (21) 

<C exp(-7 L (x,h)), (22) 

where C = L{\y\ - 1)/4. In the derivation above, (O follows since there are \y\{\y\ - l)/2 
distinct terms with s^s. For (1221) . we have used the fact that Q(x) < |exp(— ^-). 

We would like to note the similarity of (f2T|) and ([22]) to the conventional error rate expressions 
for single user wireless communication systems. Actually, the term 7 L (x, h) can be interpreted 
as a network signal-to-noise ratio (NSNR) measure that characterizes the overall performance 
of the network. 
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B. Diversity Analysis of the Relay Selection Scheme 

For K = L = 1, we have shown in [fTTTl that a feedback scheme based on relay selection can 
achieve diversity (R,0). Here, we generalize this result to any L. 

For K > 1, due to both multiuser interference and its manifestation in Theorem [fl it is not 
clear whether diversity (R, —R) would be achievable. The main goal of this section is to show 
that it is indeed achievable with a GQ that maximizes the NSNR, and surprisingly, again using 
a simple relay selection codebook. 

The relay selection codebook can be defined as C s = {e r : r = 1, . . . , R}, where e rq = 1 for 
q = r, and e rq = for q ^ r. Then, for any K and L, We define our fGQ as 

GQ Cs (h) = arg max 7 L (e r , h) , (23) 

e r 6C s 

where, for any relay selection vector e f , we have from (fl~2l) that 

2 



7 L (e r , h) = — min min 



J2 k =i( s k-h)^/PVji 



kr 



\n l 2 P 
\9r£\ -rfi 



4 -s^s- -f 1 + E f =l \f kr \2p Sk + M 2p Rr 



(24) 



Note that GQ Cs chooses the relay selection vector that maximizes the NSNR. 

In the following theorem, we show that, for both relay-broadcast and relay-interference net- 
works, GQ Cs achieves maximal diversity by finding an upper bound on the NER: 

Theorem 2. There are constants < C3, C4 < 00 that are independent of P such that for all 
P sufficiently large, 

NER(GQ Cs ) < C 3 ^, K = l, 

NER(GQ Cs ) < C 4 -|^, K>1. 

Proof: Please see Appendix ■ 

In other words, the relay selection scheme with an fGQ achieves maximal diversity for any 

K. It is remarkable that full diversity is achieved regardless of the number of transmitters and 

receivers. 

Note that our selection scheme requires |~log 2 -R] feedback bits. With [log 2 -Rol feedback bits, 

where R G {1, . . . , R — 1}, diversity orders (Ro, 0) and (R ,—R ) are achievable for K = 1, 

and K > 1, respectively, simply by considering the selection scheme for any fixed R of the 

relays and disregarding the others. 

In practical networks, we may not have a GQ that knows the entire CSI of the network. In 

such situations, we would like to characterize the achievable performance using LQ encoders 

that know only a part of the CSI. 
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V. Diversity with LQs 

In the previous section, we showed that a GQ using relay selection can achieve full diversity. 
Motivated by this result, we expect that a relay selection based LQ will achieve high diversity 
orders. In this section, we design two such LQs: An fLQ that achieves maximal first-order 
diversity, and a vLQ that achieves maximal diversity. Both quantizers will have similar structures. 
We construct them using the localization method 11201 , in which we synthesize an LQ out of 
an existing GQ. The synthesized LQ and the GQ share the same codebook. For our particular 
quantization scheme, we use the GQ GQc s in (|23T) as the basis of our LQs. Since GQc s is based 
on relay selection, all of our LQs will be based on relay selection as weljj. 

A. Localization 

Let LQi f }/ denote a generic localization of GQ Cs . For the synthesized quantizer LC)l f }/, the 
superscript indicates whether it is fixed-length (f ) or variable-length(v); and £, iV are design 
parameters that we shall specify later on. For a particular channel state h, the components of 
the synthesized quantizer operate as follows: 

1) LQ Encoders: For notational convenience, oj r £ = 7^(e r , h). The £th LQ encoder calculates 
u ri , r — 1, . . . , R. In other words, it calculates its own contribution to the NSNR for all possible 
relay selection vectors. Then, it quantizes each of the possible contributions using a scalar 
quantizer 

N(x) = J ». 3» 6 {0, . . . , W - 2} such that x 6 K, („ + 1)0, x e R (26) 

1 N, otherwise. 

Its output message is the concatenation of R sub-messages J\f(oj r i), r = 1, . . . , R. 

2) An Illustration of the LQ Encoders: Let us now illustrate the operation of the LQ encoders 
with a simple example with R = 3, and L = 2, as shown in Fig. HI For some fixed channel 
variances, power constraints, and channel state h/, suppose that u>n = 1.7, u>2i = 0.8, W31 = 1.2, 
uj 12 = 0.28, 0J22 = 0.67, and W32 = 2.3. In the figure, each of these local NSNR values are 
represented by a disk (•) on the real axis. Since we are using an LQ, u rl , r = 1,2,3 can be 
calculated only by the first receiver, and similarly, co r 2, r = 1, 2, 3 can be calculated only by the 
second receiver. Note that the GQ has access to all the local SNRS and in this example, selects 
the relay with index argmax re {i 2 ,3} miii£Ui r £ = 3. 

After the LQ encoder calculates its local NSNR values, it quantizes them using a scalar 
quantizer J\f that is uniquely determined by the parameters £ and N. In our example, we use 

N = 5 bins and set £ = |. Each bin is represented by a half open interval (E )) on the real 

axis. The output message of the LQ encoder is the concatenation of its quantized local NSNR 
values (submessages), shown as frames with a dashed outline, on the right hand side of the 
figure. 

5 ln principle, the localization method itself is applicable to any GQ with any codebook; it is not limited to relay selection 
based GQs. However, for a general GQ, it is very difficult to analytically determine the performance of the synthesized LQ. 
Therefore, we focus only on the localization of relay selection based GQs. 



15 



RXi w 2 i 

RX 2 ^22 

W32 



0.28 



0.8 



0.67 



■*■ 



Bin 1 



1.2 



4- 



Bin 2 



1.7 



* 



Bin 3 



00 



00 



2.3 



AT(w 2 i) 
AT(w 3 i) 

A/"(wi 2 ) 



■f 



Bin 4 



with r Bin 
AT = 5. I 1 f 2 
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3) Compressors: In general, there are i? sub-messages, each with iV possible values. There- 
fore, for a fixed-length synthesis LQ| n , at each channel state, each receiver feeds back |\Rlog 2 N~\ 
bits without any compression. 

For a variable-length synthesis LQ| N , we use a lossless compressor that produces an empty 
codeword (of length 0) whenever J\f (f2 rf ) = N, Vr, and otherwise a codeword of length [log 2 (A^' R - 
1)] bits that can uniquely represent each J\f(fl r{ ). In other words, for a given channel state, the 
number of feedback bits produced by any receiver is either bits or \\og 2 (N R — 1)] bitsj. 

After all the L feedback messages of the receivers are exchanged between the receivers and 
the relays, each of them decodes the feedback bits using the local decoder. The decoder operation 
is the same for each receiver and relay. 

4) Decompressor: First, a decompressor perfectly recovers all the submessages from all the 
receivers, J\f(oo r e), r — 1, . . . , R, £ — 1, . . . , L. All of these submessages are passed to the LQ 
decoder. 

5) An Illustration of the LQ Decoder: For clarity of exposition, let us first present the LQ 
decoder for the example scenario in Section IV-A2I and the same channel state h'. A more formal 
description of the general LQ decoder operation will be presented afterwards. 

If the empty codeword is not allowed, one can use a "0" (a codeword of length 1 bit) instead of the empty codeword, and 
append a "1" to each remaining codeword of length \log 2 (N R — 1)] bits. The resulting codewords are uniquely decodable as 
well. Then, all of the results in this paper will hold for the case where the empty codeword is forbidden, given that the required 
feedback rates are increased by 1 bit. 

Also, note that one can achieve a better compression by using entropy encoders instead of the suboptimal compressors that 
we employ. Even though the localization method was introduced originally with entropy encoders, the compressors that we use 
in this paper will be good enough for our purposes. 
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In general, the main goal of the LQ decoder is to imitate the GQ as good as possible. For 
our particular example, the GQ selects the relay with index argmax re { 1 2,3} u> r , where co r = 
min{a;, rl ,c<; r 2}. Then, the first goal of the LQ decoder should be to determine to r . However, the 
LQ decoder only knows the quantized local NSNR values, J\f(co r £), r = 1,2,3, £ = 1,2, as 
shown in Fig. 0] Therefore, it cannot determine the exact value of co r . However, as we shall 
describe in what follows, it can perfectly determine a subset of R where u r resides. 

For any u e R, H{oj) = n => u G [f , ^), n = 0, . . . , 3, and Af{u) = 4 => u G 
[2, tx)). We can use these facts to determine the possible locations of the local NSNR values, as 
represented in Fig. [5] by half-open intervals ( lump) of K. 
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Fig. 5: Possible locations of the local NSNRs according to the LQ Decoder. 



Since U\ = mmjwii, coi 2 }, and we know for sure that w u G [|, 2) and W12 G [0, |), we should 
have o;i G [0, |). Using the same arguments for all r, we can obtain oj 2 G [|, 1), and W3 G [1, |). 
We have thus determined the possible locations of u r , as shown in Fig. [6l by having access only 
to the quantized versions of u r . 

The LQ decoder's main goal was to find argmax re { li2j 3} u r . Using the possible locations of 
ijj r that we have found, it is now clear that the third relay should provide the best NSNR. The 
LQ decoder's output will be e 3 . Note that this is the same output as the GQ output. Therefore, 
for this particular channel state, the LQ operates in the same manner as the GQ. 

However, the LQ decoder will not be this lucky in general. As an example, another channel 
state might result in uj\ G [0, \) and ^2,^3 G [1, §). In this case, the LQ decoder will know for 
sure that both the second relay and the third relay provides a larger NSNR than the first relay. 
On the other hand, it cannot determine which one of the second and the third relays provides the 
best NSNR. Therefore, it chooses one of them, and its decision may not be the optimal one that 
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Fig. 6: Possible locations of the NSNRs according to the LQ Decoder. 



would instead be provided by the GQ. We shall quantify the effect of such suboptimal decisions 
later on. 

6) LQ Decoder: We now give the general and formal description of the LQ decoder. 

Let lZ g = {q : min^ccv = max r min^ oj T e] denote the set of indices from which our GQ in 
(1231) produces its outputjj In other words, lZ g is the set of indices of relays that provide the 
maximal NSNR. Also, let IZi = {q : J\f(mm£cu q £) = A/"(max r min^cu^)}. Note that lZ g C IZi. 
Moreover, due to the structure of J\f, not only 



Af (mm u q ij = mmAf(u q i), 



(27) 



but also 



M (max min u r p) 

v r e ' 



max M (min uj r i) (28) 

r v £ ' 

= maxminA/^Coirf). (29) 

r t 

max r mhi£.A/Xcow)}, an ^ T^i can be easily calculated by 



Therefore, IZi = {q : mmiAf(uj q £) 
the LQ decoder. 

Since lZ g C TZi, the LQ decoder can determine which relay selection vector(s) can possibly 
provide the maximal NSNR. In general, it can choose any one of the relay selection vectors that 
are indicated by TZe. But, to be more precise, we define 

LQ^ (h) = arg max e ,. gCs min^. N(u r e) ■ (30) 

7) Localization Distortion: Let us now study two possible cases of interest regarding the LQ 
output: If lZ g = IZi, then the LQ output provides the same NSNR as the GQ output. Otherwise, 
the LQ might make a suboptimal decision. This results in what we call the localization distortion 
(LD), given by 



LD(£, AT) 4 NER(LQ^ } ) - NER(GQ Cs ). 

lZ g is not necessarily a singleton, but our definition of the arg max guarantees that the GQ output is unique. 



(3D 
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A useful upper bound on the LD can be calculated as: 

l(f|v 



NER(LQ^ } ) = E h [CNER(GQ Cg (h),h)|7^ = TZ g ]?(TZi = U 6 



E h [CNER(GQ Cs (h),h)|^ ^ Tl g ]P(7li =f K g ) (32) 

< NER(GQ Cs ) + E h [CNER(GQ Cs (h),h)|7^ ^ TZ g ]P(TZi ^ K g ) (33) 

= NER(GQ Ca ) + E h [CNER(GQ Cs (h),h) | \Ki\ > \K g \\P(\Ki\ > \K g \) (34) 

<NER(GQ Cs ) + LD u (£,iV), (35) 

where LD U (£, N) is the upper bound on the localization distortion, given by 

LD U (£, AO 4 E h [CNER(GQ Cs (h),h)||7^| > 2]p(|7^| > 2). (36) 

B. Maximal First-Order Diversity with an fLQ 

Our main result concerning the fLQs is given by the following theorem: 

Theorem 3. Let £ f = log^ P, and N± = 2. Then, for P sufficiently large, the NER with LQl Ni , 
which uses a fixed R feedback bits per receiver per channel state, is upper bounded by 

NER(LQ* f>JVi ) < C 5 ^|^, K = l, 

hK 2R P (37) 

NER(LQ| fiJVf ) < C 6 ^^, K>\. 



where < C5, Cq < 00 are constants that are independent of P. 

Proof: Please see Appendix [Ql ■ 

In other words, using a fixed R feedback bits per receiver per channel state, we can achieve 
diversity (R, -R) for K = 1, and diversity (R, -2R) for K > 1. Since (R, -R) < (R,0) for 
the broadcast network, and (R, —2R) < (R, —R) for the interference network, our fLQ has a 
second-order diversity loss compared to the optimal performance for both types of networks. 
Also, it is straightforward to show that, using R bits, where R G {1, . . . , R}, we can achieve 
diversity gains (i? , — Ro) and (R , —2R ) in relay-broadcast networks and relay-interference 
networks, respectively. 

The scalar quantizer resolution for our fLQ is log 2 N f = 1 bit per local NSNR. In what follows, 
we show that, by appropriately increasing the resolution with P, one can achieve maximal 
diversity, while the compressors make sure that the feedback rate remains bounded. 

C. Maximal Diversity with a vLQ 

For vLQs equipped with entropy coding, we have the following result: 



19 



Theorem 4. Let e > be a fixed constant that is independent of P. For any A that satisfies 



< e < A < P, let £ v - x 

AT 



\, and 

[AlogA + PAlogP + 1], 



A log A + PA log 



P 



logP 



K = l, 
K > 1. 



(38) 



TTzen, for P sufficiently large, we have 



LD u (£ v ,iV v )<CV 
LD u (£ v ,iV v )<C 8 



1 



AP R ' 
log R P 



K 



K > 1, 



(39) 



an J, zn addition, the feedback rate of the 



R K L Qg v ,Af v ) < Cio 



receiver satisfies 
IokP 



P ' 
log 2 P 
_ P~ ' 



K 



K>1, 



(40) 



where < Cj,Cs,Cg, Ciq < oo are constants that are independent of A and P. 

Proof: Please see Appendix E ■ 

We now describe several consequences of this theorem for K > 1. The consequences for 

K = 1 will be analogous. 

Let us first recall from ([35]) that NER(LQ£ >JV ) < NER(GQ Cs ) + LD U (£ V , N v ). We have found an 

upper bound for NER(GQ Cs ) in Theorem[2] An upper bound for LD U (£ V , N v ) is given by Theorem 

IH Combining the two bounds, we have NER(LQ| >N ) < (C 4 + CsA" 1 ) 1 -^^- . In other words, our 

vLQ achieves maximal diversity. 



Moreover, using the same arguments as in the previous paragraph, we have NER(LQ| 



N) 



< 



NER(GQ Cs ) + ^ l0 pu P ■ Thus, by increasing A, the array gain performance of our vLQ can be 
made arbitrarily close to the one provided by the GQ, at any finite power level P. 

What is more interesting is the behavior of the upper bound on the feedback rate given by 
(l40l) . As P grows to infinity, the required feedback rate decays to zero. In other words, both the 
diversity and array gain benefits of NER(GQc s ) can be achieved using arbitrarily low feedback 
rates, when P is sufficiently large. 

VI. Simulation Results 

In this section, we present numerical evidence that verifies our analytical results. We assume 
that each receiver attempts to decode all the symbols from all the transmitters. In other words, 
V e = {l,...,K}, W. In the graphs, "GQ" represents GQ Cs in (J23]), "fLQ" denotes LQ| ;JVf with 
£j and Nf as defined in the statement of Theorem |3l Also, "vLQ-A" represents LQ? Nv that is 
uniquely determined by the parameter A as in the statement of Theorem HI 
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A. Networks With Equal Parameters 

In Fig. [7J we show the performance results for a network with K = R = L = 2, er? = 
a <L = Pn r = Ps k = 1, Vr, k, I, and <Si = S 2 = {+1,-1}. For this network, the NERs with the 
GQ, fLQ, and vLQs for A = 2~ 15 , 2" 12 , . . . , 2 12 , 2 15 is presented in Fig |7al The horizontal and 
the vertical axes represent P in decibels (dBs), and the NER, respectively. 

We can observe that both our GQ and vLQs achieve the maximal diversity (2, —2), while the 
fLQ achieves diversity (2, —4). Moreover, as we increase A, the array gain performance of our 
vLQs can be made arbitrarily close to that of the GQ. 

In Fig. |7bl we show the SERs at the first receiver for the same network. The horizontal axis 
represents P in decibels, while the vertical axis represents the SER at the first (second) receiver. 
As a result of our choice of network parameters, the SERs of each receiver is the same. Also, 
a particular quantizer achieves the same diversity as in Fig. [Taj On the other hand, since the 
SER is upper bounded by the NER, any quantizer in Fig. |7b] provides more array gain than it 
does in Fig. [Taj Indeed, due to the symmetry of the network parameters, the SER performance 
is around 1.6dB better than the NER performance for all quantizers. 

The corresponding feedback rates of our vLQs are shown in Fig. [7c] The horizontal axis 
represents P in decibels, while the vertical axis represents the feedback rate of the first (second) 
receiver in bits per channel state. Similarly, due to our choice of the network parameters, the 
feedback rates of each receiver will be the same. We can observe the validity of Theorem HI 
as for any A, the required feedback rate decays to zero at high P. Also, by increasing A, the 
performance of the LQs can be made arbitrarily close to the one provided by the GQ, while still 
using very low feedback rates. As an example, at an NER of 1CT 5 , vLQ-2 15 needs 1.25 bits per 
channel state per receiver on average and performs only 0.25dB worse than the GQ. At a SER 
of 1(T 5 , vLQ-2 6 uses 0.65 bits, and GQ performs only 0.8dB better. 



B. Networks With Unequal Parameters 

Our results also hold for networks with unequal power constraints and/or channel variances. 
To demonstrate that, we consider a network with K = R = 3 and L = 4. The parameters 
of the network are assumed to be p Sl = 1, Ps 2 — 1-3, Ps- A — 0.7, p Rl = 0.6, p R2 = 2, 



PRs 



>2 



0.7, Si = S 3 = {+1,-1}, and S 2 = {e j ^ e : 6 e {!,..., 4}}. Also, we assume that 



a 



J k r 



Fkr, V. 
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Fig. 7: Performance results for a network with K = R = L = 2. 
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(c) Feedback rates. 

Fig. 7: Performance results for a network with K = R = L = 2 (continued). 



In Fig. M we show the NERs with the GQ, fLQ, and vLQs for A = 2~ 15 , 2~ 12 , . . . , 2°, 2 3 . 
The results are analogous to what we have observed in Fig. [7a] Both the GQ and the vLQs 
achieve the maximal diversity (3, —3), while the fLQ achieves diversity (3, —6). Moreover, as 
we increase A, the array gain performance of our vLQs can be made arbitrarily close to that of 
the GQ. 

The SERs at the first and the third receiver are shown in Fig. [8b] and Fig. [8c] We can observe 
that, unlike the previous network with equal parameters, the SERs at each receiver is different 
for this network with unequal parameters. In particular, Fig. [8b] reveals rather counterintuitive 
results: The fLQ outperforms the GQ at low P, and some of the vLQs provide a higher array 
gain than the GQ. The reason of these behaviors is that the GQ is optimized with respect to 
the NER, which takes into account the SERs of all the receivers. Therefore, as far as the SER 
at a receiver is concerned, one cannot claim the optimality of the GQ. For the NER, the GQ 
outperforms all the other quantizers, as shown in Fig. [8a] 

For the vLQs, the feedback rates of the first and the third receivers are shown in Fig. [8d] 
and Fig. [8eJ For both figures, the feedback rates decay to zero as P grows to infinity, verifying 
Theorem [4] 
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(b) SERs at the first receiver. 



(c) SERs at the third receiver. 
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(d) Feedback rates at the first receiver. 
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Fig. 8: Performance results for a network with K = R = 3, L — 4. 
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VII. Conclusions and Discussions 

We have studied quantized beamforming in wireless relay-interference networks with any 
number of transmitters, receivers and amplify-and-forward (AF) relays. Our goal has been to 
minimize the probability that at least one user incorrectly decodes its desired symbol(s). 

We have introduced a generalized diversity measure in order to have a more precise description 
of the asymptotic performance of the network. It has encapsulated the conventional measure as 
the first-order diversity. Additionally, it has taken into account the second-order diversity, which 
is concerned with the transmitter power dependent logarithmic terms that appear in the error 
rate expression. 

First, we have shown that, regardless of the quantizer and the amount of feedback that is used, 
interference results in a second-order diversity loss in our network model. Care should be taken 
though when making a general statement, as in this work, we have focused on AF networks 
with a short-term power constraint. For other forwarding methods, such as decode-and-forward, 
the diversity results may be different. Even under the restriction of using AF relays, one can 
use a long-term power constraint and achieve higher diversity. Also, the side information at the 
relays may be exploited for a better performance, though we believe this will not improve the 
diversity. 

Second, we have designed a relay-selection based global quantizer (GQ) that can achieve 
maximal diversity. Then, using our GQ and the localization method, we have synthesized fixed- 
length and variable-length local quantizers (fLQs and vLQs). Our fLQ has achieved maximal 
first-order diversity. Our vLQ has provided not only maximal diversity gain, but also an array 
gain performance that can be made arbitrarily close to the one provided by the GQ. Moreover, it 
has achieved all of its promised gains using arbitrarily low feedback rates, when the transmitter 
powers are sufficiently large. 

Regarding the LQs, there are many open problems that we have not addressed in this paper. 
One important problem is to determine whether there exists an fLQ that can achieve maximal 
diversity. Another goal might be to generalize our relay-selection based localization result to 
show that any GQ can be localized to synthesize an LQ that can achieve the same array gain as 
the GQ. Due to the complicated nature of our distortion functions, the latter goal seems difficult 
to accomplish, even though we have observed its validity by simulations. 

Appendix A 
Upper Bounds on the PDF and CDF of Vt r 

First, let us present some useful lemmas. 

Lemma 1. Let f c and f, be zero-mean real Gaussian K x 1 random vectors, with equal diagonal 
covariance matrices E[f c fJ] = E[f s fJ] = K, Ka > 0,Wi, Kij = 0, Vi ^ j, and zero cross- 
covariance E[f c fJ] = 0. Let f = f c + jf s denote the complex Gaussian random vector with real 
and imaginary parts given by f c and f s . Also, let X = |(s, f)| 2 , where s G C K — {0} is a fixed 



25 



vector, and W = ||f || 2 — X. Then, there is a constant < A < oo, such that for all x > and 
w > 0, we have 

\ K 

fx,w(x, w) < ° ^~ 2 exp [-\ (x + w)] . (43) 

Proof: Let f = Uf for a unitary matrix U that satisfies e x = Us. Also, let X { = \fi\ 2 , and 

X=[X 1 ---X K ] T . Note that 

X = X 1 = |(s,f)| 2 = |(e 1 ,f)| 2 , (44) 

and since U is norm-preserving, 

W = ||f || 2 - X = ||f || 2 -X l = E?=2 X *- (45) 

Now, let Y\ = X\, Y 2 = W = J2i=2 -^-i, and Y k = X k , k — 3, . . . , K. Using such a transforma- 
tion of RVs [23 1, we have 

V2 yi-yz V2-T, k J 3 Vk 

/n,K 2 (2/1,2/2) = / / ••• / /xd/i,!/2-Ef=3fc?/3, ■■■,yK)dy 3 ---dy K . (46) 
00 

In the following, we find an upper bound for /x(x) for any x = [x% ■ ■ • xk] t with Xi > 0, Vi. 

Let U c = 9ft(U), and XJ S = S(U). The real and imaginary parts of f can be calculated to be 
f c = 9ft(f ) = U c f c — U s f s , and f, = S(f ) = U c f s + U s f c . Then, it is straightforward to show that 

K cc 4 E[f c fJ] = U c KUj + U s KUf , (47) 

K ss ^ E[f s fJ] = K cc , (48) 

and 

K cs ^ E[f c fJ] = U c KUf - U s KUj. (49) 

Therefore, K cc and K ss are symmetric matrices, and K cs = — Kj s . The latter implies that for 
any x e M^, x T K cs x = 0. Using these facts, we now show that K cc + jK cs is positive definite. 
For any x e C K , we have 

x"(K cc + jK cs )x = (xj - jxJ)(K cc + j'K cs )(x c + jx s ) (50) 

= X c xv cc X c X c xv cs X s + X s xv cs X c + X s xv cc X s + 

j(xjK cs x c + xjK cc x s - xf K cc x c + xf K cs x s ) (51) 

= xf K cc x c + 2xf K cs x c + xf K cc x s (52) 

= xf E[f s fJ]x c + 2xf E[f c f J]x c + xf E[f c fJ]x s (53) 

= E[(xff s + xff c ) 2 ] ^ ^ (54) 

= E[((xJU s + xJU c )f c + (x^U c - xf U s )f s ) 2 ] (55) 

= xf Kxi + x£Kx 2 , (56) 
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where xi = Ujx c + U^x s , and x 2 = Ujx c — Ujx s . But, 

Ux = (Uj + jUf )(x c + jx s ) = x 2 + j Xl 



(57) 



and since Ux ^ 0, either xi ^ or x 2 7^ 0. Also, since K is positive definite, either xj Kxi > 0, 
or x 2 r Kx 2 > 0. Thus, x H (K cc +jK cs )x > 0, Vx G C K — {0}, and K cc +jK cs is positive definite. 
Let A + jB = (K cc + jK cs y\ and % = tan _1 (Sij/A^). According to flM Eq. (24)], the 
joint PDF of Xj, i = 1, . . . , K can be expressed as 



/x(x) = (4vr)- K det(A + jB 



exp 



;/(*) 



d</>, 



(58) 



where 



K 



A' 



/(#) = ^ Ai^ + 2 ^ {A% + Bj)* v^i cos (^ - & + %) 



»j/' 



(59) 



i=l 



i<j 



and D = A + jB is a Hermitian matrix flM Eq. (21)]. 

Since /(0) is continuous, and the range of integration [— n,n] K is a compact subspace of 
R K , there exists </>* e M^ with cf>* = [<j>\ ■ ■ ■ <j>* K ], such that f (</>*) < /(</>), V0 G [-vr, n] K . As 
a result, 



jx(x)<2- K det(D)exp - ~f(<f>* 



Now, let 



X c 
X s 



xa-cos(^)] t , 
i^sin(0^)] T . 



'Xl cos^-j ■ 
/xi" sin(0*) • • 

Then, using (l59l , /(<£*) can be expressed as 

/(**)= K(xLDx C0S + xJ in Dx 



(60) 

(61) 
(62) 

(63) 



We have shown that K cc + jK cs is positive definite. It follows that D = (K cc + jK cs )~ l 
is also positive definite, and thus has eigenvalues Aj > 0, VI Also, since D is a Hermitian 
matrix, it admits a spectral decomposition D = J2i=i ^j u « u f > where Uj, i = 1, . . . , K form an 
orthonormal basis for C A . It follows that 



Ji 



X cos-'-^ X cos — X C0S JJX C0S 



> Allx 



E i= i Ai(uf x 
2 



cos/ 



where A = minj Aj. Similarly, we have 



x Dx ■ > Allx ■ 

-^■sin sin -^ /v |l sin| 



(64) 
(65) 

(66) 



27 



Using (|65l) , (|66l ) and (|63l ), a lower bound on /(</>*) is given by 

/(0*)>K(A||x cos || 2 + A||x sin || 2 ) (67) 

= A £?!*<, (68) 

Then, using (|68l) and (1601) , we can find an upper bound on /x( x ) as 

/x(x) < 2- A 'det(D) Uli exp(-f^) (69) 

^AfniiexpC-AoxO, (70) 

where A = §. For the last inequality, we have used the fact that det(D) = Yli=i A» — A K • 
The lemma follows by substituting (|70l) to (|46|) and performing the integration. ■ 

Lemma 2. Lef X 1; . . . , X n &e n non-negative possibly dependent RVs, and Z = min n X n . Then, 

n 

fz(z)<J2f^), (VI) 



and 



F z (z)<J2F Xi (z). (72) 



»=i 



Proof: Let us recall Leibniz's integral rule: For functions of a single variable a(z), b(z), 
and of two variables f(x,z), we have 

b(z) b(z) 

§~ z I /(*, ^ = / f dx + /(&(*), ,)g - /(a(*), *)^. (73) 

0(2) a(z) 

Note that CZ2]) easily follows from (f7fl . We thus first prove (fTTb . Let Z fc = miii{I b . . . , X fe }. 
We will show that fz k {z) < Y^,i=ifxi{z), f° r an Y 1 < A; < ra by induction. For fc = 1, it 
is obvious. Suppose it is true for n > k > 1. We have fz k (z) — 12i=i fxX z ) ■ Noting that 

-Zfc+i = min(Z/ c ,X/ c+ i), 

fz h+1 (z) = fz k (z) + fx k+1 (z) - -^F Zk ,x k+1 [z, z) (74) 

< E /*(*) - ^ J J h k ,x k+1 (n,v)dudv (75) 

4=1 

= ^fxAz) - I < ^ / /z fc ,x fc+1 (M,w)dt; \ du ~ fz h ,x k+1 (z,v)dv (76) 

4=1 o I o Jo 

= 5Z -fa (*) ~ / fzk,x k+1 (u, z)du - I f Zk ,x k+1 (z, v)dv (77) 

4=1 

fc+1 
<£/*(«), (78) 



i=i 
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where for both (1761 ) and (1771) , we have used Leibniz's integral rule. This proves (I7T1) . Integrating 
both sides of (1711) from to z proves (1721) . ■ 

We can now find the desired upper bounds on the PDF and CDF of Q r . 

Proposition 2. For all ui > 0, y > and P sufficiently large, 

1) IfK = l, 

faM < Cnexp^-Cia^VoCw), (79) 

and 

i^fo + u) - Fa r {y) < CnwVo(w), (80) 

4M = -^(l+^), (81) 

anJ < Cn, C12 < oo are constants. Otherwise, 

2) lfK>\, 

/n» < C 13 exp (-C14J) VH, (82) 

F n „(y + w) - i^fo) < C 13 fwV(") + ^Z/ 2 ) , (83) 

and in particular, for y — 0, 

Fn r (w) < Ci3W^(w), (84) 

vv/zere 

logP / 1 1 1 1\ 

^(w) = ^^(l + u; 2iogp+a/ logP p), (85) 

anJ < C13, C14 < 00 are constants. 

Proof: First we prove the case for iT > 1. Let £l r ,t,s,a — 7£ S s( e »->h)- Note that Q, T = 
min^ iS ^ s tt r ,i„ Sl §. First, let us first find an upper bound on the PDF and CDF of Q r ,i,s,s- 

Consider a fixed r, I, and s/s. For notational convenience, let us define Z = n r ^ )S) g. From 
(124)) . we have 

1 Sfe=l( S fe - h)\/P~S~ k fkr 



\n l 2 P 



(86) 



Now, let us rewrite (l86l) in a more compact form. First, we define 

? = lVP^flr---VP^fKr} T , (87) 

S±[si-s 1 ---s K -s K ] H , (88) 

X^|(f',<5>| 2 , (89) 

V = \g r e\ 2 p Rr , (90) 

W = ||f|| 2 -X, (91) 



29 

where P = f '/ j| f || . Then, we have 

aXYP 2 

v _ (92) 

" l + XP + WP + YP' v ' 

where a = \\\S\\ 2 . Using a transformation of RVs [23J, the PDF of Z can be expressed as 

oo oo 

f(\ f ff f z(l + xP + wP) \ axP 2 (l + xP + wP) f 

Hz) = J J fv [ axP ,_ zP ) {axP 2_ zPr f*,w(x, w)dxdw. (93) 



Now, let A y = (cTg re p Rr )~ l . Substituting the PDF of Y, and using Lemma [Q we have 

oo oo 

,/w V* f f ( A y z(l + xP + wP)\ 
h{z)< *™ / /exp - V _„ 9 — M x 



r(ir-i) 77 *a axP 2 - z p 






aaP 2 (l + xP + wP) _ A ^_ 2 -a «,, ,^n 

. ^ Q — - — -e A ™ x dxw K 2 e A ™ w dw, (94) 

[axP 1 — zP) 1 

where oo > A xw > is a constant that is independent of w, x, and P. The inner inte- 
gral can be evaluated first by a change of variables u = axP — z and then using the facts 
that r o °° x v - 1 e-^ x -^(dx = 2(p/j)*K v (2y/p T j),P,'y > fl26l Eq. 3.417.9], and K^(x) = 
Ki(x), Vx G K 11251 Eq. 9.6.6], respectively. Then, after some straightforward manipulations, we 
can rewrite (l94l) as 

oo 



T(K-l) J V u ; a 2 P 2 V 7e A ™» 



where k = ^jA\^ N \ y z(a + z + awP)/(a 2 P 2 ), and A = (A xw + \ y )/a. It follows that 

_ oo 

fz(z) < ^K-'l)^ /fa" 1 "^) + g-Vifo(/c))t^- a e-^«dti;. (96) 

o 
Now, let us find an upper bound for K (k) in (|96| ). According to [|25l Eq. 9.6.24], we have 
K„(z) = f^° e~ zcosht cosh(ist)dt, t,v6l Moreover, since cosh(z/£) is an increasing function 
of z/, ./f„(z) is also an increasing function of za It follows that 

K (k) < K v (k), u > 0. (97) 

Also, from ifTTl Eq. 25], we have 

KJk) < 2 u - 1 T(u)k- u , u>0. (98) 



Now let us set < v < 1. In this case, 



T(u) = v' x Y{y + 1) = v- 1 / e~* dt < is' 1 / e _t dt = v~ x . (99) 

o o 
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Combining (1971) , (1981) and (1991 ) gives us the desired upper bound 

#o(«) < 2 1/ - 1 z/ _1 «- I/ . (100) 

Using (fTOOl) and the fact that ^(k) < 1 031 Eq. 25], ([96]) can be further bounded as 



/z(*) < 



A^T 1 exp(-A- 

r(K-i) 



(XP- l + 2 



1 i ov-l..-l„-l. 2-i/\ „..K-2-\ xw w, 



-l„-A-§r 



< d 5 ^ e~ A ^ 



'// : 'ir J )w K ~' z e~ x ™ w <lw 



(101) 
(102) 



P -i + p-2+^^-5 u 1 + wP + z y-2 W K - 2 e - x ™ w dw 



where C 15 = A*" 1 ^ - l)]" 1 max{A, 2A xw A y max{±, £}, 2(A xw A y max{i, £})£}. Also 
since (x + y)" < x^ + y", Vx, y G M, < v < 1, we have 



fz(z) < C l5 v- l e- x ? 



,1-£d1-£ , ~l-Z\„..K-2-\ xw w, 



P- L + p- 2+u z-z / (1 + w 1 ^? 1 - + ^-2 )w K - 2 e- A ™ w dw 





(103) 



Ci 5 e" 



wP 



>-i+*vi , „i-S 



i + z - 5P -^(i + ^-*)r(jr - i)A xw ^ +1 + 



-I (K-0A^ + M 



(104) 



< 2C 15 r(K) max{l, A~*+\ A^K^^P- 1 ^ (i + z -\ + z^P' 1 ) (105) 

= Ci 6 -^e- X * (l + z~***y + ^^TSipp- 1 ) (106) 

where Ci 6 = 2eCi 5 T(K) max{l, A~^ +1 , A~^}, and we have substituted v = -^-p to obtain 
(TT06T) . 

In general, the constants C*i6 and A in (11061) depend on r, £, s, and s. Let Ci6,r,^s,s and A r ,£, Si s 
denote the dependent versions of Ci 6 and A, respectively. Using Lemma [2] we have 



w 



< X^ eX P ("CM-) (l + *"^ + Z^^P- 1 



(107) 



(108) 



where Ci 3 = 2Za5^|(| J^l — 1) max rAs ^s Ci 6iri i Si s, and Cu = min^^g A r ^ jSj s. This implies the 



upper bound on the PDF of fi r in the statement of the lemma. Finally, using (1 1071) , 

f Qr (u)du 
y 



< 



< 



< 



CislogP 



4 


P 


C13 


logP 


2 


p 


C*13 


logP 



/ I 1 + X^^r + X 1_I ^P — ] dx 



u + (y + lu) 21 °s 



P 



a; + u 21 °g p + 



P _ y 1 - 21ogP _|_ 

y 2 ~\^p 4- 2w 2_I °g : P 



(!/ + w) 



1 
p 



1 
p 



<c'j, ( w^(w) + -^V >°^ 
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(109) 
(110) 

(HI) 

(112) 
(113) 
(114) 



where we have used Holder's inequality and the fact that (y+co) a < y a +oj a , y,z>0, 0<a<l 
for (111 31 ) . This concludes the proof for K > 1. 



For K = 1, let X r = |/i r | 2 Psi» Y r = p Rr min £ |^| 2 , and T r 



i+x r p+y r P' 



From (124]) . we 



have Vt r = CyjT r , where Cyj = |min Sl ^ Sl (si — §i). 

Now, note that X r ~ r(l,p 5l a 2 lr ), and F r ~ r(l,p Br (E < 0" 1 )- Therefore > X, ^n r = 
1, . . . , R are independent exponential RVs with finite variances. The PDF of T r with such X r 
and Y r is given by ifTTl Eq. 22]. Using [11, Eq. 28] without omitting the exponential function, 
and noting that /h r (u;) = ^-fr r (^-), we can show that (1791) holds. 

Finally, (f8~0l) follows (up to a constant multiplier) from (1791 and (11091) . This concludes the 
proof. ■ 



Appendix B 
Proof of Theorem [I] 

We start with a lower bound on the CNER. By definition, we have CNER(x, h) > SERj ML (x, h) , V 
Suppose that, for some k G T>i, a genie reveals all the transmitted symbols but s k to the £th 
receiver. The error rate of this genie-aided scheme provides a lower bound on the CNER. Without 
loss of generality assume that 1 G T>\, and let us fix some s l5 s 2 G «Si with s\ ^ s 2 . We have 

CNER(x,h) > ^Q( v /27 u (x,h)), where 

ER p I — T I 9 I ' / |9 n 

r=l hrVPr9rlX r \ \Sl ~ S 2 \ Z Ps 1 



7 U (x, h) 



4(l + Ef=lPrbriri' 



2U |2"| 



(115) 
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Let us find an upper bound on 7 u (x, h) for any h, and x 6 X. We have 



7 n (x>h) <l^-^)>/^E r=1 /ir^^ r | 2 (n6) 



4(l + EtiP;l^i 



2\ x \2\ 



< I Si ~ g2| 2 -Psi I Er=l hry/fi.y/R9rlX r \ 2 
4 i2 + Etl^brl| 2 |Xr| 2 



I c o I 2 P 



Si 






4 Eti(i + p'Ai| 2 )|5vl 2 ' (118) 

where x = Ax is the projection of the beamforming vector onto the hypersphere with norm 
R. Applying the Cauchy-Schwarz inequality to (11181) . and then using the fact that p' r < p r , we 
have 

\A A |2p R I f |2_ D|„ 12 

7°(x,h> < if- - ; z ' p » x: ^ f^f; aw 

If K = 1, we use the following upper bound that follows from (111 9b . 



+ p r R\g rl \ 
|s'i -s 2 | 2 P 5l 



7 P (x,h)< |ai ^ rSl J2\f lr \ 2 . (120) 

r=l 

This upper bound is, up to a constant multiplier, the same as the SNR of a maximal ratio 
combining system with R branches. The error rate of such systems is known to be lower bounded 
by a constant times P~ R , as stated in the theorem. This concludes the proof for K = 1. 
For K > 1, we use (11191 ) to further bound 7 u (x, h) as 

R 2 \6 6 I 2 If \ 2 \n l 2 P P 

7 x, h) < max - -—- (121) 

4 r ! + l^k \fkrrPs k + R\g r irPRr 

P 2 lc e I 2 If |2 1 „ |2p p 

< ^ |gl ~ 521 max l/irl fri ggA 122 

4 r l + l/irl^x + l/iJrl^Sa+^rll 2 ^ 

< fi 2 |ii-g2| 2 m a x r {4 r < lPgl p Rr } X r F r P 2 

" 4min{l, a^ m ,cr^ & , ^p^} r 1 + X r P + W r P + Y r P ' 

where X r . = vj 2 f lr , Y r = v~ 2 g r i, and W 7 ,. = a7 2 f 2r - Note that X r , F r , W r ~ r(l, 1) and they 
are independent. Let Cis denote the constant multiplier in (11231) . and Z" = (X r F r P 2 )/(l + 
X r P + W r P + Y r P). Thus, we can rewrite (fT23l) as 7 u (x, h) < C 8 max r . Z r u . Now, let 

Z u = max r ^, (124) 

A l 



NER L (Q) £ 1 ^ f E[Q(v^COT U )]. (125) 

Since NER(Q) > NER L (Q), VQ, it is sufficient to find a lower bound on NER L (Q). Using the fact 
that Q(x) > ^Y^e 2 , we have 



NER L (Q) > — \-= [ ^^exp(-zC 18 )f z v(z)dz. (126) 

\Sl\y/1T J l + ZZ 
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We thus need to find a lower bound for the PDF of Z u . Using order statistics, we have 

R R 

fz<z) = Y,fz»M)X{Fz l {z). 



(127) 



r=l 



9=1 



In the following, we find a lower bound on the PDF and CDF of Z", for any r. We first evaluate 
the PDF of Z®. Using a transformation of RVs |f23l , it can be expressed as 



fz?(z) 




,. z(l + xP + wP)\xP 2 (l + xP + wP) _ x , 
fy r I — — — = — - ) ^^ — — -e dxe dw. 



o o 



X p 2 -zp J (xp 2 - zpy 



(128) 



This PDF is in the same form as (l92l) in Proposition [2l and can be evaluated using the same 
methods discussed therein. We have 



oo 



U « 1+ z p+z W J « l+ Z p+ ' ) » 



2(1 + 2z + iwP) ^ /4z(l + w? + z) 



P 2 



P 2 



e-"dw 



(129) 
(130) 



>exp(-|)^/ 2 (l + W P + ^ ( ^ 1+ ; 2 P + ^ )M w) d W . (131) 
o 

Using the fact that for any z > 0, X"„(z) = - log(f) ~7e + (1 - Te)^ + (1 + 1/2 - l e % + ■ ■ ■ 
11261 . we have K (z) > — log(|) — 7 e , and thus 



oo 



log 



z(\ + z + wP) 

P 2 



27e 



e~ w dw (132) 



{oo 
(-2 7e + 2 log P - log z) f{l + z + wP)e- w dw- 

oo -. 

{l + z + wP) log(l + z + wP)e- u, dw I 

exp( - ^)^{(-2 7e + 21ogP - logz)(l + z + P)- 
1 + z) log(l + z) + Plog(l + z) + P + Pe^E/ ' 



(133) 



(134) 



Using the facts that logz < log(l + z) < z, and 



e^Ei ( — 1£ ) < log ( 1 + ^— ) < log(l - ^ + P) < log(l -h z + P + zP) 



P 



1 + z 
log(l + z) + log(l + P) < 2 + log(2P) < z + 1 + logP, (135) 
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2z\ 1 



we can show that 

/^(*0>exp(-^)-^{piogP + 2bgP(l + *)- (136) 

[(2 7e + z){\ + z + P) + z(\ + z) + 2(z + 1)P] } (137) 

> exp( - ^-IjpiogP - P[(2 + z){2 + z) + z(l + z) + 2(* + 1)] } (138) 



2z\ 1 



log P - (2^ 2 + 7,2 + 6) 



After some straightforward manipulations, (11391 ) leads to a more compact lower bound 



where 



Mz) > <f>(z) 



e-^-[logP-14(l + z 2 )]. 



2x\ ] 

I— [logP- 14(1 + x 2 )] dx 



For the CDF of Z v r , we have 

z 

Fz*{z) ^ / ex p( 



> exp (- ^)|; [log P- 14(1 + ,2 2 )] 

> z6(z). 



PJP 

2z\ z 



(139) 

(140) 
(141) 

(142) 

(143) 
(144) 



We can now find a lower bound for the PDF of Z u . Suppose that P > e 14 , and let z — 
("~u~ ~~ 1 ) 1 - Then, 0(z) > for z < z , and <fi(z) < 0, otherwise. Using (11271) . for z < z , it 
follows that 



Mz) > Rz 



H-lj.fi/ 



= P^" 1 exp {—) p\ [ lo S P - 14 ( X + *' 
>Pz^ 1 exp(-2P^)— s [log P- 14(1 + z 



2M« 



2M R 



pi? 



P^ 1 exp(-2P^)-^ ^ rM log^ 4 P (-14)* (1 + 

i=o ^ z ' 



2\i 



> P^- 1 exp(-2P2)^[log iJ P - f^ i L> ^ ' - ■-"-' 



i=l 



\ og R - 1 PW 1 (1 + z 



2\i 



> Rz 11 ' 1 exp(-2Rz)— w [\og R P - R2 2R ~ l U R log^ 1 P(l + z 



2R 



pR 



(145) 
(146) 

(147) 
(148) 

(149) 
(150) 



Since fz a {z) is a PDF, fz^{z) > 0, \fz. Therefore, for z > z , we can choose any negative 
function as a lower bound on f&>(z). But, (11501 ) is negative for z > z . Thus, it is a lower bound 
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on fz u (z) that holds for all z. We can therefore use it to bound (11261) as 

oo . 

L/ , R log R P f z R ~2 7ic .„«, 

BEE (Q) '- iW^J TT2? ' " ,dz " 



oo , 

- % , i 4 ° g 5 / - ( +Z V ^+^dz. (151) 

o 
The first integral in (11511) can be lower bounded by 

oo 1 1 1 

/" _f^ e -*(Ci 8 +2fl) d2 , > f ^^ e - z ^ +2 ^dz (152) 

i l + 2 2 -J l + 2z 





e -(C 18 +2R) 



> / z R -^dz (153) 



o 

2 e -(Ci8+2i?) 



3(2P+1) 
For the second integral in (11511) . we have 



oo , 

v fi-i/i i ^2R. 







r(p+§) r(3P + § 



(154) 



l Z ^1 + i ) e - 2 (C 18+ 2R) dz < />"!(! + JS 2fl )e -*(C i8 +2fl) d? (155) 



- + ^-T- ( 156 ) 

(C 18 + 2i?)^l (C 18 + 2Rf R +\ 
Substituting (fT54l) and (fT56l) to (fT5Tl) . it follows that 

NER L (Q) > 2C 19 p- /? (log i? P-C 20 log i? - 1 P), (157) 

for some constants oo > C w , C 20 > independent of P. 
Finally, C 20 log 72 " 1 P < \ log R P, VP > exp(2C 20 ), and thus 

NER L (Q) > C 19 P- R log fi P, (158) 

for all P > exp(max{14, 2C 20 }). This concludes the proof. ■ 

Appendix C 
Proof of Theorem [2] 

We provide a proof for K > 1. The proof for if = 1 is very similar. Thus, we skip it for 
brevity. 
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Let f2 = max r il r , where Vl r = 7 L (e r ,h), as defined in Appendix lAl Then, we have 

NER(Q Cs ) < C E[exp(-fi)] 



/ri JX 

r=l <7=1 



vAU))] 1 ' X do; 






00 

<RC C& f u R ~ l e- w (l + 



U 21ogP + UJ - logP J ^UJ 



R 



.R- 



/v 1 



< m^CoC^ J u^e* ( 1 + u'^^ + u/ 1- ^-^- I <L- 



log R P 



r(p) + r p- 



/? 



21ogP 



+ r 2i2 - 



p 



\orPJ P R 



<c 



21" 



(159) 
(160) 

(161) 

(162) 
(163) 

p R > ( 164 ) 

where C 2i = P3 R-1 C C^|max{r(2P),r(|)}, and (1 1601) follows from the order statistics of 
independent RVs. For (11611) and (|162l) . we have used Proposition [2l and Holder's inequality, 
respectively. This concludes the proof. ■ 

Appendix D 
Proof of Theorem [3] 

Let us prove the theorem for K > 1. The proof for K = 1 is very similar. It is thus omitted. 

Let n and Cl r , be as defined in Appendix O We need to find an upper bound on the localization 
distortion. According to (l36l) . it is sufficient to calculate the CNER given \1Zi\ > 2. Note that 
\1Zi\ > 2 if and only if there exists r,q E {1, . . . ,P}, r ^ q such that Af(Cl r ) = J\f(£l q ) = 
jV(fi). Depending on Af(fl), we divide the calculation of LD u (£,iV) to two separate parts as 
LD u (6,iV) = EL 1 LD«(e,iV). 

The first part is concerned with the case J\f(Q.) = 0, or equivalently, tt r E [0,£ f ), Vr. Since 
the decoder chooses one of the P relay selection vectors, the NSNR is at least min r il r . Using 
Proposition [2l we have 



V SI 

LD"(£ f , 2) < C / • - • / exp (- minojr) Y[ fa r (u r ) J} du} r 

r 



2R 



< C (PC; 



13 J 



M n P 



pR 
r 2_R 



l + (PlogP) 2l °sP + (PlogP) 1 lo « F — 



< C (PC 



13 



r WP ( 1 



pR 



+ e^ + 



PlogP 
P 



(165) 

(166) 
(167) 

(168) 
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\og 2R P 
< C 2 2^pj^, (169) 

for a constant < C 22 < oo, and all P sufficiently large. 

For the second part, we consider the case M{Vt) = 1 -<=>- 3r G TZ, fi r G [£ f , oo). In this 
case, the minimum NSNR is £ f , and we simply have LD"(^f , 2) < C -^n. 

Combining the final upper bounds for the two parts, we have LD u (£ f , 2) < C23 og pR for some 
constant < C23 < 00. This concludes the proof. ■ 

Appendix E 
Proof of Theorem [4] 

We prove the theorem for K > 1. The proof for K = 1 is very similar and skipped for brevity. 
Let f2 and Q r , be as defined in Appendix O Also, for simplicity of notations, let £ = £ v , 

iV = A v , and 5 = (iV v - l)f . Depending on Af(Q), we divide the calculation of LD U (£, A) to 
three separate parts as LD U (£, A) = ^=1 l&Kt* N )- 

The first part is concerned with the case where A/"(f2) = •<=>■ Q r G [0, £), Vr. In this case, 
the NSNR is at least mm r f2 r . Using Proposition |2l we have 

LD?(£,A) <C f ■■■ /exp (-nunw r ) JJ/OrCwr)!!^ ^ ^^ V(0- 0™) 


Now we consider the term ^>(£) in (11701) . For future reference, we shall calculate an upper bound 
for the more general quantity given by ip(n£), for any n G {1, . . . , N — 2}. We have 

logP 



^«) 



P 



1 + (n£) 21 °g p + (n^) 1 lQ s p — 



(171) 



WP / !_ _L_ 1 L_ 1 1 \ 

= — ^— ll + n 2i og p\2io S p +n i logp^iogp j (172) 

log P / 1 1 1 \ 

< — |- I 1 + A^p 4. n^Aisip j (173) 

<!^(i + vs +en? i), (1 74) 

where the last inequality follows from A < P. Moreover, for all n G {1, . . . , N — 1}, n£ < 
(A - l)f < log A + PlogP - P log log P + i < e- 1 + (P + 1) log P. Combining with (fl74l) . 
we can argue that there is a constant < C 24 < 00 such that 

V(n£) < cJ^A (175) 

for all P sufficiently large. Using (fTTOl) . it follows that LD?(£, A) < C(,( fjf 24 lo f,* F . 

For the second part, we evaluate the cases for which J\f(Q) G {1, . . . , A — 2}. For each 
n G {1, . . . , A — 2}, suppose that i > 2 of fi r are in the interval [n£, (n + 1)£), and the rest 



3S 



R — i of them are in [0,n£). The minimum NSNR is at least n£. Also, there are (^) possible 
ways to choose which Q r will be in [n£, (n + 1)£). Therefore, 

Ar _ 2 < n£ (n+1)? (n+l)£ 



a n-2 7 " ? ^y^ ^"T, 

ld$«,ao<£ j; j] /.../ / ... / 

i=2 fCejr R n=l n n c t 

;^t^-j U U n£ n£ 



_R-i integrals j integrals 

R 

exp(-mino; r ) JJ/n r (^ r ) J^[ dw r ]^[ day, (176) 



r=l 



rS/C r'g/C c 



where J^ is the collection of all possible z-combinations of the set {1, ...,R} (e.g. J^ 2 3 = 
{{1,2}, {1,3}, {2,3}}), and /C c = {1, . . . , R} — /C. Then, similarly, we can use Proposition [2] to 
arrive at 



LD^,iV) 

R N-2 

*= 2 /CeJK^ n=1 



i?,-i 



mO + ^H) 2 -^ 



p 2 



t=2 



/? 



JV-2 



< CoC^E ( ■ P ~* E nfl ~ ie ~^ fl_i ( n ^ 



n=l 



eV(0 + ^(i + ^ 2i ) 



log' P 

p2i 



U<r R P - 



P R 



i=2 



7V-2 



i?+i 



7V-2 



(£ R + £ R - { P- 1 ) Y^ n^e-^ + 5__ Y n R+i e~ ni 



n=l 



n=l 



(177) 

(178) 
(179) 



where (11781) follows from Holder's inequality, and the fact that (n^) 2l ^ T °sP < (l + n 24 £ 2 *), i > 1. 
For (11791 ), we have applied (11751) . Now, we shall evaluate the summations with respect to n in 
(11791) . The following lemma provides a useful upper bound: 



Lemma 3. Let f be a non-negative real valued Riemann integrable function with f(x) < 
oo, Va; e K. that is increasing on the interval (— oo, b), and decreasing on (6, oo). Then 

(180) 

n=n\ 



ri2 £ 

J2f(n)< J f(x)dx + 2b. 



hi 



Proof: Let n b = [b\ be the largest integer less than b. Assume that n\ < n b < n 2 . Then 



n b -l n b -l n+1 n b -\ "+ 1 J 

E /(») = E / f^ dx - E / f^ dx = / • 

n=ni n=ni n=ni " 

x * n L n n\ 



f(x)dx, 



(181) 



where the inequality follows from the fact that / is increasing in the range of integration. Also, 



"2 



"2 



n+1 



"2 



n+1 



i^2 /» "<2 /» /• 

E f^= E / f( n ) dx < E / /(«-l)dx= / f(x)dx, (182) 

n=n b +2 n=n 6 +2 n=n b +2 



n b +l 



V-) 



where the inequality follows since / is decreasing on [n b + 1, oo), and thus for n b + 2 < n < 
x < n + 1, f(n) < f(x - 1). Finally, combining (11811) and (11821) . 

n 2 n b — l n 2 

n=ni n=«i ra=nj,+2 

rij, — 1 ri2 

< J f (x)ix + b + b+ J f {x)dx (,84, 

«1 "b + 1 

II) 

<Jf (x)ix + 2 b , (,85, 

"1 

which is the desired inequality for ni < n b < n 2 . The other cases can be proved similarly. We 
skip them for brevity. ■ 

Note that the function f(x) = x l exp(— x£) has a global maximum at x = i/£ with f(x) = 
i l £~ l exp(-i). Moreover, for any < a < b < oo, f f(x)dx < / °° f(x)dx = T(i + l)^^ l+1 \ 
Using Lemma |3j for any i < R, we have 

JV-2 

j2 n R-i e -< < T (R - i + i)^-(«-i+i) + 2(i? - $)*-*£-(«-*) (i 86 ) 



71=1 



< r(i? - i + i)^-(»-*+i) + 2(i? - j)*-*£-(*-i+i) (187) 

<2( J R-i) /? - i r (i?_m) - (188) 

where the second inequality follows from the assumption that £ < 1. 
Using (11881) . (11791 ) can be bounded as: 

i R p -R 
LD^, iV) < C (2C' 1 3) ii -^- ^ [2(P - i) R - i ^- 1 +r 1 P- i ) + 2(R + i^C'P'*] (189) 

i=2 

< C (2C 13 ) R ^^ E W^ " 2)^ 2 r ' + 2(2i?) 2R f- 1 ] (190) 

j=2 

< PC (2C 13 )* [4(i? - 2) R - 2 + 2(2P) 2R ] ^^ (191) 

where the second inequality follows from the assumption that P > £ -1 . 

For the last part, we consider the cases for which N(Z) — N — 1. The minimum NSNR is 
(N- l)f = S > log A + P bg(j^p), and we have LD|(£,iV) < C e^ < C ^#. 

Combining the final upper bounds for LDy(£, N), i = 1, 2, 3, LD?(£, JV) < C 25 ^j£ for all P 
sufficiently large, and a constant < C25 < 00 that is independent of P and A. This proves the 
upper bound on the LD. 
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Finally, by the definition of our compressor in Section IV-A31 we have 

R,(LQ^J = |Tog 2 (JV B - l)]P(3r, Q r£ < N) (192) 

R 

< \R\og 2 N]J2 p (ttri<Z) (193) 

r=l 

logP 



<RC 2 Al + R\og 2 



AlogA + i?Alog( ) +2 

loeP 



P 



(194) 



log 2 P 



< C 26 ^^, (195) 

for some constant < C 2 % < oo, and P sufficiently large. The third inequality follows from 
(11751) . This concludes the proof. ■ 
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